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ABSTRACT 
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1^-^ Aims. We investigate the structure of the Galactic magnetic field in the 2nd Galactic quadrant using radio continuum polarization data from the 
f*"*) I 325 MHz WENSS ( WEsterbork A'orthern 5ky Survey) survey in combination with earlier single-dish observations. 

^—i Methods. We study gradients in polarization angle along Galactic longitude and latitude in the region 130° < / < 173° and -5° ^ b < 
\^ 31°. These gradients were determined with a new method that we developed to efficiently and reliably fit linear gradients to periodic data 
like polarization angles. Since the WENSS data were obtained with a synthesis array they suffer from the 'missing short spacing' problem. 
We have tried to repair this by adding an estimate of the large-scale emission based on the single-dish (Dwingeloo) data obtained by Brouw 
, and Spoelstra. Combining all available data we derive a rotation measure map of the area, from which we estimate all 3 components of the 
I magnetic field vector. 

O Results. In the part of WENSS where large-scale structure in polarized intensity is relatively unimportant, we find that the magnetic field is 
predominantly perpendicular to the line-of-sight, and parallel to the Galactic plane. The magnetic field components along the line-of-sight and 
^ along Galactic latitude have comparable values, and the strength of these components is much smaller than the strength of the total magnetic 
. . , field. Our observations also cover part of the so-called 'fan' region, an area of strong polarized intensity, where large-scale structure is missing 
^ ■ from our WENSS data. We tentatively show that Faraday rotation occurring in front of the Perseus arm is causing both the WENSS RM and 
the RM towards the fan region observed in previous single-dish surveys, and we suggest that the fan is formed by local emission that originates 
in front of the emission we see in WENSS. 

Key words. Magnetic fields - Radio continuum: ISM - Techniques: polarimetric - ISM: magnetic fields - ISM: structure - Techniques: 
miscellaneous 



1. Introduction 

The diffuse radio emission from our Galaxy provides a pow- 
erful tool for studying the various components of the interstel- 
lar medium. For example, the well-known all-sky survey by 
Haslam et al. ( 198 1 ) at 408 MHz was used by Beuermann et al. 
J1985> to construct a two-component model of the synchrotron 
emissivity of our Galaxy. Large parts of the Milky Way have 
been studied also at other frequencies both in total intensity 
and in polarized intensity, e.g. Berkhuijsen J1971> . Reich et al. 
SW l I, Uyaniker et al. (2003 1, Bingham & Shakeshaft l fT967t . 
Brouw & Spoelstra ( 1 976 ) . 

High-resolution polarization studies came about with the 
work by Junkes et al. ('1987'), and an aperture synthesis ar- 
ray was first used by Wieringa et al. Cl993j . Haverkorn et al. 
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(I2000I I were to the authors' knowledge the first to use an inter- 
ferometer to derive information on the magnetic field structure 
in detail using the diffuse emission of our Galaxy. Presently, the 
International Galactic Plane Survey (IGPS) project is study- 
ing a major part of the Galactic plane up to a couple of de- 
grees on either side and at several frequencies (the polarization 
data from the IGPS are discussed in Taylor et al. 2003 for the 
Canadian Galactic Plane Survey and by Haverkorn et al. l2006al 
for the Southern Galactic Plane Survey). 

Equipartition arguments have been used to estimate the 
strength of the magnetic field (Beck et al. l2003t . To study the 
topology of the field various groups are using different tech- 
niques. In particular pulsars and extragalactic sources have pro- 
vided a large amount of information on the large-scale field and 
on possible reversals in the direction of the field (see e.g. Han 
et al. 2004>.2006i. and Brown et al. 2003 ) . Extragalactic sources 
sample the entire line-of-sight through the Galaxy, whereas 
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pulsars can tell us about variations along the line-of-sight if 
we know their position in the Galaxy. Most pulsars are how- 
ever confined to the Galactic plane, a drawback extragalactic 
sources do not have. 

One major problem with both extragalactic sources and 
pulsars is that except in parts of the Galactic plane (Brown et 
al. .2003. Haverkorn et al. .2006bJ the surface density of extra- 
galactic sources and pulsars is low, leaving large gaps between 
the sampled lines-of-sight. As a result only large-scale struc- 
ture (LSS) in the magnetic field can be inferred. Extragalactic 
sources can also have an intrinsic rotation measure, which 
makes it more difficult to extract the Galactic contribution from 
the observed rotation measure. 

Pulsars have previously been used to provide information 
on the small-scale magnetic field (Rand & Kulkarni 1989 Han 
et al. I2004> . The polarized diffuse Galactic radio background 
can also be used for this purpose, as shown by Haverkorn et 
al. (I2004> . who derived the strength of the magnetic field com- 
ponents along the line-of-sight and perpendicular to it, mak- 
ing it possible to study both the small-scale and the large-scale 
magnetic field. One problem with this method is that the dif- 
fuse emission is not very strong, and it can be severely depolar- 
ized along the line-of-sight. Correcting for this is complicated 
by the fact that the amount of emission and depolarization can 
vary along the line-of-sight. Furthermore the diff'use emission 
fills the entire telescope beam, which makes modelling more 
involved than for pulsars and extragalactic sources which are 
bright point sources. 

In this paper we study the large-scale properties of the 
polarized diff'use Galactic background at 325 MHz using the 
WENSS survey (Rengelink et al. '1997V We selected a region 
spanning about 1000 square degrees between 130° < / < 173° 
and -5° < b < 31°. This area encompasses the Auriga and 
Horologium regions studied by Haverkorn et al. ( I2003all2003bl 
using multi-frequency data, as well as a large part of the so- 
called 'fan' region, an area bright in polarized intensity, where 
the magnetic field component in the plane of the sky is very 
regularly distributed (Brouw & Spoelstra il976. and Spoelstra 
[19841 . 

In Sect.|2we describe the WENSS dataset that we use and 
the reduction of this dataset. Since this dataset was obtained 
with an interferometer, LSS will be missing from our observa- 
tions. In Sect. |3] we estimate the contribution of LSS, and we 
add this estimate to WENSS. We present the results from the 
original WENSS dataset and that including the LSS estimate in 
Sect.0] In Sect.|5]we describe how we analyse the polarization- 
angle data. We also discuss the robustness of our LSS recon- 
struction. In Sect. |6l we present the results of our analysis in 
terms of large-scale gradients along Galactic longitude and lat- 
itude. In Sect.EJwe compare the gradients in polarization angle 
that we obtain in areas where no LSS is missing from the in- 
terferometer observations to gradients in the rotation measures 
derived by Spoelstra ( 1984i. This enables us to estimate how 
reliably gradients in polarization angle can be translated into 
gradients in RM. In Sect. |8]we discuss the physical picture we 
can draw from our observations; in particular we attempt to re- 
construct all components of the 3D magnetic field vector. In 



this section we also consider some of the implications of our 
observations for the nature of the the fan region. 

2. Description of the data 

2.1. WENSS 

WENSS, the Westerbork Northern Sky Survey, is a low- 
frequency radio survey of extragalactic sources at 325 MHz 
above S = 30° that also contains a wealth of diff'use polariza- 
tion data (Rengelink et al. HWTb . WENSS has a 5 MHz band- 
width divided over 7 channels, and its sensitivity is 1 8 mJy (5cr) 
at 325 MHz. The survey was carried out using the Westerbork 
Synthesis Radio Telescope (WSRT). The WSRT consists of 14 
25m dishes, 10 of these are fixed and 4 can be moved along a 
track, which improves (u,v) coverage. 

To map large parts of the sky in a reasonable amount of time 
WENSS uses the mosaicking technique in which the telescope 
cycles through 78 or 80 pointings (depending of the declina- 
tion of the mosaic), integrating each pointing for 20 seconds, 
and then moving on to the next. If the pointing centers are on 
a square grid with a separation of 1.5° the off-axis instrumen- 
tal polarization is very effectively suppressed (down to the 1 % 
level - Wieringa et al. 1993 1. In this way each mosaic covers 
roughly 200 square degrees, and in every night each pointing 
can be observed about 12 times, yielding visibility data along 
12 'spokes' in the (u,v) plane. 

Each WENSS mosaic is formed by combining six 12hr ob- 
serving runs. In each run the telescopes sample the (u,v) plane 
with a 72m baseline increment, and different runs have a dif- 
ferent shortest baseline ranging from 36m to 96m in 12m inter- 
vals. This results in baselines from 36m to 2760m. At 325 MHz 
this gives a maximum resolution of 54" x 54" esc 6 (FWHM). 
In the analysis of the data we used a Gaussian taper with a value 
of 0.25 at a baseline length of 250m to increase the signal-to- 
noise ratio for more extended structures. This degrades the res- 
olution to 6.7' X 6.7' CSC S. By adding six 12hr observing runs 
the first grating ring is at 4.4°, i.e. outside the 3°x3° area that 
we mapped in the individual pointings from which we construct 
the mosaics. 

The 6 mosaics we selected were mostly observed at night. 
Observing at night limits solar interference, and strong iono- 
spheric RM variations are considerably smaller when the ob- 
servations are not taking place during twilight. Information on 
the mosaics can be found in Tabled For the observations that 
started in the afternoon ionospheric effects are probably small 
since polarized point sources are still point-like in the data. 
After the polarization calibration described in the next section 
we regridded the Stokes Q and U maps from the equatorial to 
the Galactic coordinate system using the AIPS task REGRD, 
after which we calculated the polarized intensity and polar- 
ization angle. Polarization angles were coiTected for the (lo- 
cal) parallactic angle between the equatorial and Galactic north 
poles. 
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Table 1. Summary of the observations. For each mosaic the dates (yy/mm/dd) and starting times (UT) are shown of the individual 
12 hour observing runs that make up the mosaic. Each run is indicated by the length of its shortest baseline. For each mosaic we 
also show the central coordinates (both equatorial and Galactic), as well as the resolution (FWHM) after applying a Gaussian 
taper. We used this resolution in further analysis of the data. 



Mosaic 


WN50.074 


WN50.090 


WN66_045 


WN66_064 


WN66_083 


WN66_102 


(,a, ();ic)5o.o 


/ 4 , J\j ) 


[yyj , J\J ) 


{^■J , 00 ) 


[OH , 00 ) 


{OJ , 00 ) 


,00 ) 


(m 


(157°,6°) 


(163M4°) 


(136°,8°) 


(142M2°) 


(146M9°) 


(149° ,26°) 


Resolution 


6.9' X 9.0' 


6.9' X 9.0' 


6.9' X 7.6' 


6.9' X 7.6' 


6.9' X 7.6' 


6.9' X 7.6' 


36 m 


92/01/07 


92/01/13 


93/10/18 


92/01/12 


92/01/09 


92/01/08 




15:24 


16:06 


18:45 


14:26 


17:36 


17:14 


48 m 


91/12/31 


92/01/06 


93/10/24 


92/01/05 


92/01/02 


92/01/01 




15:52 


16:33 


18:22 


14:53 


16:21 


17:41 


60 m 


91/12/24 


91/12/30 


93/10/28 


93/01/10 


91/12/26 


91/12/25 




16:20 


17:01 


18:06 


14:31 


16:49 


18:09 


72 m 


91/12/06 


91/12/02 


93/11/01 


91/12/01 


91/12/05 


91/12/04 




17:41 


18:51 


17:50 


17:11 


18:12 


19:32 


84 m 


91/12/10 


91/12/16 


93/11/07 


91/12/15 


91/12/21 


91/12/11 




17:48 


17:56 


17:27 


16:16 


17:44 


19:04 


96 m 


91/12/17 


91/12/23 


94/10/09 


91/12/22 


91/12/19 


91/12/18 




16:47 


17:28 


19:22 


15:48 


17:16 


18:37 



2.2. Polarization calibration 

Rengelink et al. ( 1997 1 carried out the total intensity calibration 
as part of the WENSS survey. For details about this part of the 
data reduction we refer the reader to their paper Polarization 
calibration requires additional reduction steps that we carried 
out using the NEWSTAR data reduction package, and we will 
briefly describe these here. For a proper (mathematical) treat- 
ment of the different steps involved in the polarization calibra- 
tion we refer the reader to the articles by Hamaker et al. (.1996 1 
and Sault et al. ([1996). 

Corrections for deviations in the alignment and the elliptic- 
ity of the antenna dipoles, which cause total intensity 1 to leak 
into Stokes U and V, were found using the unpolarized calibra- 
tor sources 3C48 and 3C147. The flux scales of these sources 
is set by the calibrated flux of 3C286 (26.93 Jy at 325 MHz), 
determined by Baars et al. (I1977> . 

Due to an a-priori unknown phase offset between the re- 
sponses of the X and Y dipoles there is cross-talk between 
Stokes U and V. In principle this phase difference can be cal- 
ibrated with a polarized calibrator source, but since such a 
source was not available for the WENSS observations we deter- 
mined the correction by assuming that Stokes V contains only 
noise, and minimizing the amount of signal in V. For the 5 mo- 
saics observed in January 1992 the correction we found in this 
way was -9°, a demonstration of the stability of the WSRT on 
timescales of weeks, whereas for the WN66_045 mosaic, which 
was observed almost 2 years later, it was -f26°. 

Finally we have to correct for the different amounts of iono- 
spheric Faraday rotation during the different nights in each mo- 
saic and between mosaics. In each mosaic we found 2-3 bright 
intrinsically polarized extragalactic sources, and we used the 
variations in the polarization angle of these sources over the 
different nights to align the polarization vectors. We did not 
correct for variations in ionospheric Faraday rotation within 



each individual 12hr observing run; we only took the average 
polarization angle over each night. To correct for polarization- 
angle offsets between mosaics we compared polarization an- 
gles in the area of overlap between mosaics. The corrections 
we found gave a closure error of about 1° when going from 
WN50.090 to WN50_074,WN66.064, WN66.083 and back to 
WN50.090. The WN66.045 and WN66.102 mosaics only have 
one neighbouring mosaic, which means that closure errors can- 
not be investigated for these mosaics. We assume that the align- 
ment of these mosaics is of the same quality, even though the 
width of the distribution of polarization-angle differences be- 
tween the WN66_045 and WN66_064 mosaics is almost a fac- 
tor of 2 larger compared to other mosaics. 

3. Estimating large-scale polarization structure 
missing from WENSS 

If an interferometer does not provide visibilities for the central 
part in the (u,v) plane LSS in the sky will be missing from the 
observations. For the WENSS data structure larger than about 
1.5° - 2° is missing. We estimated the importance of structure 
on these scales from single-dish data from Brouw & Spoelstra 
(1976 'BS') and Spoelstra (1987"). This (linear) polarization 
dataset was obtained with the Dwingeloo 25m telescope at 
frequencies of 1411, 820, 610, 465 and 408 MHz, beamsizes 
(FWHM) range from 0.6° at 1411 MHz to 2.3° at 408 MHz. 
These data were however sampled on an irregular grid, with 
an average pointing distance of about 2°-3° below a Galactic 
latitude of 20°, and about 5° above a latitude of 20°. 

Before we can combine the BS and WENSS datasets we 
have to resample the BS data onto the WENSS grid, and we 
have to extrapolate the BS data from 408 MHz to the WENSS 
observing frequency of 325 MHz. These issues will be ad- 
dressed in the next two subsections. In the third subsection we 
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discuss how we added the two datasets. In Sect. l5.2l we will dis- 
cuss the robustness of the maps we made of the WENSS data 
including our estimates of the missing LSS. 

3.1. Interpolating the Brouw & Spoelstra 
measurements 

We addressed the problem of the irregular BS sampling by con- 
volving the BS datapoints with a Gaussian kernel. The severe 
undersampling of the BS data sets a lower limit to the kernel 
width, and a reasonable upper limit can be found by requir- 
ing that the interpolated data should not be too smeared out. 
We used 5 different kernels, with FWHM ranging from 3.1° to 
8.7°, to interpolate the BS measurements of Stokes Q and U 
and polarization percentages at their 5 observing frequencies. 
In addition we interpolated their maps of the intrinsic polar- 
ization angle of the emitted radiation and the RM maps they 
derived. 

3.2. Extrapolating from 408 MHz to 325 MHz 

The second problem is how to get from the BS measurements, 
with a lowest observing frequency of 408 MHz, to the 325 MHz 
of the WENSS data. To predict the polarized intensity at 325 
MHz, Ppied,325, we start with the observed polarized intensity 
at 408 MHz, PobsAoa, convert it to total intensity and scale it up 
to 325 MHz using a power-law, and finally convert this back to 
polarized intensity at 325 MHz: 

Ppi-ed325 = PchsA)& X (408/325)'" X ppt,,, (1) 

where ppf^c is defined as 

_ polarization % at 325 MHz 
" polarization % at 408 MHz 

For the polarization angle at 325 MHz, 0pred,325, we calculate 
the amount of Faraday rotation from: 

<^pred,325 = 0O,BS + ^■'^BS'i^325MHz) (3) 

where /1(325mhz) is the observing wavelength corresponding to 
325 MHz. The (pp^ed are wrapped back to [-90°,90°]. For these 
predictions we use the intrinsic polarization angle of the emit- 
ted radiation 0o and RM as derived in Spoelstra 1 19841 ('(/in rs' 
and 'RMbs'), together with the power-law index they adopt for 
the brightness temperature, 2.7 (which gives the 0.7 in Eqn.^ 
since the 2.7 is used for scaling brightness temperatures, and 
we convert total intensities), and the polarization fractions they 
derive using this power-law at the different frequencies. Eqn.|3l 
requires an interpolated map of polarization angles, but inter- 
polation cannot take the periodicity of the data into account. 
We tried to minimize this effect by interpolating ^o.bs data, 
which showed less ±90° transitions than other polarization- 
angle data. 

For the actual prediction of the large-scale Stokes Q and U 
at 325 MHz from the 408 MHz observations we must assume 
a polarization percentage at 325 MHz. From Fig. 6 in Brouw 
& Spoelstra ( 1976) it is clear that the average polarization per- 
centages decrease slightly between 610 MHz and 408 MHz for 



b=0°— >-i-20°. By using the same polarization percentage at 325 
MHz and 408 MHz (at each frequency one value for the en- 
tire mosaic) we can be sure that we are not underestimating 
the predicted polarized intensity of the BS data at 325 MHz or, 
equivalently, the importance of LSS missing from WENSS. 

We checked how well this scheme works by using the BS 
polarization data and interpolated ppfr^c maps at 1411, 820 and 
610 MHz to predict the observed polarized intensity at 820, 
610 and 408 MHz respectively. If we do not use a pixel-to-pixel 
correction for the variation in polarization percentage but just 
one correction factor for the entire mosaic, we predict on aver- 
age a polarized intensity that is about the same as the observed 
intensity. The spread (standard deviation) in the ratio of the pre- 
dicted vs. the observed polarized intensities (about 0.2) is very 
similar to that obtained with a pixel-to-pixel correction for vari- 
ations in polarization fraction, except at 408 MHz where it is 
0.5. The equivalent width at 408 MHz when the pixel-to-pixel 
correction is used is still about 0.3. This increased spread can 
be explained by varying amounts of Faraday rotation between 
pixels, which becomes more important at lower frequencies. 
If part of the emission is actually thermal Bremsstrahlung, we 
would be overestimating the amount of emission at lower fre- 
quencies by using a synchrotron power-law with spectral index 
2.7. Apparently this problem is not very important, given the 
results of our comparison between the predicted and observed 
polarized intensities. 

From these results we conclude that our recipe to predict 
polarized intensity works sufficiently well that we can attempt 
to predict polarization data at 325 MHz from the BS 408 MHz 
data. 

3.3. Adding the LSS estimate to WENSS 

For an excellent review of different methods to add missing 
LSS to interferometer observations see Stanimirovic ( I2002> . 
We used the 'linear combination approach', where the Stokes 
Q and U maps from WENSS are added to the reconstructions 
of Stokes Q and U from the BS data. However, as we could 
not cross-calibrate the WENSS and BS data the intensity scales 
are not guaranteed to be identical. A potentially more serious 
problem is that the polarization vectors in WENSS and BS may 
have different zero-points. In Sect. I5.2l we will investigate the 
effects of such errors. 

In the remainder of this paper we will us the term 'WENSS 
LSS' when we refer to the WENSS dataset complemented by a 
BS reconstruction of the missing large-scale Q and U structure 
at 325 MHz. 

4. Discussion of the maps 

In Fig.^we show the polarized intensity (in Kelvin) and polar- 
ization angle for WENSS and WENSS LSS respectively. Two 
of the brightest features in both WENSS and WENSS LSS, viz. 
the V-shaped feature at (l,b) - (161°, 16°) and the ring-like 
structure at (137°, 7°), were analysed in detail by Haverkorn 
et al. (2003a, "2003b I using multi-frequency observations. The 
most striking feature in Fig.^is the clear stratification of polar- 
ization angles in Galactic latitude. In polarized intensity there 
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170 160 150 140 130 170 160 150 140 130 

Galactic longitude [ " ] Galactic longitude [ " ] 

Fig. 1. Polarized intensity (top row) and polarization angle (bottom row) mosaics for WENSS (left column) and WENSS includ- 
ing our reconstruction of the missing LSS ('WENSS LSS'; right column). The beamsize of the observations (not indicated) is 
6.7' x6.7'csc 6 (FWHM). Properties of the reconstruction (see the text): width of the convolution kernel is 5.7° FWHM, pph-ac-^, 
A^-Q°. The polarized intensity scale in both mosaics is the same. The nested rings in polarization angle in WENSS LSS at 
Galactic latitudes b > 20° is an artefact due to a single large RM in the BS data surrounded by much smaller RM. 



are also linear features aligned with the Galactic plane, but 
these are not as clear as the structures in polarization angle. 

After adding the BS data to WENSS the ring-like structure 
that was already clearly visible in polarization angle now also 
becomes more conspicuously ring-like in polarized intensity. 

In WENSS LSS there is an area between 136° < / < 152° 
and -3° ^ ^ ^ 13° where the polarization angles apparently 
show no clear stratification with Galactic latitude. This region 
is part of what has previously been referred to by some as 
the 'fan' region, a large area of high polarized intensity with 
a strongly aligned projected magnetic field (see e.g. Figs. 6a to 
6e in Brouw & Spoelstra , 1 976. and Fig. 3a in Spoelstra 1984 1. 
Where the BS signal is strong (in the upper left and lower right 
of the maps) there are clear differences between WENSS and 
WENSS LSS. In some parts this structure in BS seems to erase 
the stratification of polarization angle we noted in the previ- 



ous paragraph, but in large parts of WENSS LSS the stratifica- 
tion is still very clear. We will return to this point in Sect. (S] 
The nested imprint in polarization angle in WENSS LSS above 
b = 20° is caused by a single line-of-sight in the BS data that 
has a large RM, suiTounded by much smaller RM. This creates 
steep gradients in RM which by Eqn.l^are translated into steep 
polarization angle gradients in the reconstructed BS data. There 
also appears to be a general change in scale of the polarization 
features above b = 20° in WENSS. In the WN66.102 mosaic, 
which is the main source of data above b > 20°, an unusually 
large fraction of the (u,v) datapoints at the shortest baselines 
had to be flagged. This could be partly responsible for the ob- 
served change in structure. However, flagging of these data in 
other mosaics did not significantly alter the maps. 
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interval] 

bin i 




best estimate of slope 



Fig. 2. Overview of how we analyse our polarization data. First 
we derive the average polarization angle in each bin (polarized 
intensity weighted), (pun- Bins are aligned with the 'x' and 'y' 
Galactic coordinates. We then fit a gradient to all (p^n in an 
interval, taking into account the periodic nature of the polariza- 
tion angles, and we repeat this for all intervals in a strip along 
Galactic longitude or latitude. The gradients can then be used 
to investigate the average gradient per strip, as shown in the 
third row. 

5. Quantifying tlie large-scale stratification of 
polarization angle: technique 

The striking stratification of polarization angle in Fig. ^ sug- 
gests that there is a clear gradient in the direction of Galactic 
latitude, with only a relatively small variation in Galactic longi- 
tude, on which small-scale structure is superimposed. In order 
to quantify this large-scale gradient it is necessary to filter out 
these small-scale modulations, and also to coiTectly handle the 
180° ambiguity of the polarization angle. 

Resolving the 180° ambiguity by minimizing the difference 
in polarization angles of consecutive datapoints is sensitive to 
small-scale structure and noise. By using longer intervals in 
longitude or latitude we can fit the large-scale behaviour in po- 
larization angle without being influenced too much by small- 
scale structure, and at the same time resolve the 180° ambiguity 
for all fitted datapoints together. 

5.1. Determination 

The different steps in our analysis of the polarization-angle data 
are illustrated in Fig.|2] To determine the large-scale gradients 
we first calculate a weighted average of the polarization an- 
gles over an area of a fixed size ('bin', the average itself is 
indicated by '0bin' in Fig- ID- For strips along Galactic latitude 
these bins are 28'x7' in size (about 4x1 WENSS beams, taking 
into account the orientation of the beam in the Galactic coor- 



dinate frame) and for strips along Galactic longitude 14'xl3' 
(or about 2x2 WENSS beams). In this way all bins contain 
about the same amount of beams, and the different shapes of 
the bins reflect the difference in steepness between gradients 
along Galactic longitude and latitude. If the spread of the po- 
larization angles within the bin is larger than 20° the bin is dis- 
carded. Furthermore if a cut is placed on e.g. the signal-to-noise 
ratio of the pixels in a bin, at least 50% of the pixels should be 
usable (not flagged), otherwise the bin will also be discarded. 

We then make a linear fit to all usable 0bin in an interval of 
a specified length, as shown in the second row of Fig. |2] We 
only made fits to intervals that contain at least 6 usable bins. 

The fits were made using a standard minimization. 
Instead of determining for every possible combination of gradi- 
ent and offset which 180° 'flips' of the datapoints minimize the 

for that combination it is also possible to determine which 
configurations of 'flips' are allowed by the data. The best fit 
to each configuration can then be found using a standard x^ 
minimization, and the configuration with the lowest x^ value 
will give the overall best fit. This method can be shown to be 
more reliable and time-efficient than an approach that 'probes' 
which gradient fits the data best. See appendix A for details. 
This technique can also be extended to situations where the 
horizontal distance between datapoints varies, as is the case 
when polarization-angle observations at different wavelengths 
are used to derive RM. 

For every strip we also fitted gradients that are shifted by 
half an interval length to prevent losing information (Nyquist 
sampling). In the remainder we include the gradients in these 
shifted intervals in our analysis. 

5.2. Robustness of the gradients including the LSS 
estimate 

When combining WENSS with our reconstruction of the large- 
scale structure there are several free parameters that needed 
to be chosen: the FWHM of the convolution kernel for the 
BS data, the possible misalignment between the WENSS and 
BS polarization vectors (A(p) and differences in the amplitude 
scales of the WENSS and the BS data (which can be ab- 
sorbed into differences in ppfmc)- We have investigated the ef- 
fect of varying these parameters by comparing histograms of 
polarization-angle gradients. In the remainder of this section 
we will consider only gradients along Galactic latitude and we 
test for only one interval length of 1.4°. As there is much more 
structure in the polarization-angle gradients along Galactic lat- 
itude as compared to gradients along Galactic longitude, gra- 
dients along Galactic latitude are more sensitive to changes 
in the parameters we are investigating. An interval length of 
1 .4° gives the largest number of usable gradients in our dataset. 
Since at least 6 bins per interval must be usable to fit a gradi- 
ent, the number of gradients fitted over short interval lengths 
will be limited. On the other hand for an area of a given size 
the number of intervals that can be fitted in that area decreases 
when the interval length goes up. 

Note that we only investigate the statistical properties of 
the polarization angles that we find after adding the large-scale 
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structure. We do not claim that individual pixels will have the 
correct polarized intensity or polarization angle. 

Taking A0=O° we found no significant differences in the 
distributions of the polarization-angle gradients for all widths 
of the convolution kernel > 5° (FWHM). Smaller widths are 
simply too narrow to smoothly interpolate the undersampled 
BS datapoints, leading to artifacts. As there are hardly any dif- 
ferences for the different widths of the interpolation kernel we 
will be using the reconstruction made using the FWHM=5.7° 
kernel for further analysis. Since the WENSS data is missing 
information on scales > 1.5° - 2° this means that we do not 
have information on scales from about 1.5° to about 5° after 
adding the BS reconstruction. 

Since the WENSS and BS datasets do not overlap in the 
(u,v) plane it is not possible to determine if there is a misalign- 
ment between the WENSS and BS polarization vectors. We 
appUed polarization-angle offsets of 30°, 60°, 90°, 120° and 
150° to the BS reconstruction with a kernel width of 5.7° be- 
fore adding it to WENSS, and found no significant change in 
the gradient distribution. We therefore keep fixed at 0°. 

Finally we considered the effects of differences in the in- 
tensity scales between WENSS and the reconstructed BS data. 
These consist of two contributions, the assumed polarization 
percentage at 325 MHz, which is important for ppfr^c (cf. Eqn. 

0, and differences in the flux scales between the WENSS and 
BS datasets (mentioned in Sect. 13. 3> . 

In our reconstruction of the LSS missing from WENSS we 
assumed that the ratio of the polarization fractions at 325 MHz 
and 408 MHz (Eqn.O was 1 . This ratio is more likely less than 

1, maybe as low as 2/3, which would imply that there is less 
LSS missing from WENSS than we assumed. In the top two 
panels of Fig. 13 we show the gradient distributions for these 
two values of p p^.^^. . The only major difference between the his- 
tograms is the fraction of gradients that are in the peak around 
radians/deg. The second peak around -2.3 radians/deg hardly 
changes. In Sect. 15.31 we identify these peaks and explain the 
difference for the two values of ppfac- For now it suffices to 
say that this difference does not influence our results. 

Since the influence of calibration differences between the 
WENSS and BS flux scales is probably smaller than the range 
of the two values of ppt^c we probed, we did not consider this 
difference in flux scale as a separate issue. 

5.3. Dependence on the value of pptmc 

In Fig. 13 we plotted the distributions of the gradients fitted to 
the WENSS LSS dataset for ppf^c^l, 2/3, 1/3 and for the orig- 
inal WENSS distribution (ppbac-O)- In this way the gradual 
change in the distribution between the different reconstructions 
becomes clear. 

Gradients in WENSS are found in two peaks at +2.3 ra- 
dians/deg and -2.3 radians/deg. Since WENSS LSS is con- 
structed as the vector sum of the WENSS and BS datasets, the 
relative strength of the BS and WENSS polarization vectors 
determines if the gradients in WENSS LSS show the clear bi- 
modal behaviour of the original WENSS dataset. In this way 
the negative mode and tentative positive mode in the distribu- 
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Fig. 3. Histograms of gradients fitted to the WENSS LSS data 
and to WENSS. The top panels show the result for a reconstruc- 
tion of WENSS LSS using a ppfac (defined in Eqn.|2ji of 1, 2/3 
and 1/3 respectively. In the bottom panel we plot the distribu- 
tion of the gradients fitted to the WENSS data without the LSS 
reconstruction. In all cases the width of the convolution kernel 
used to reconstruct the LSS is 5.7° (FWHM), A0=O°, and the 
interval length over which the gradients were fitted is 1.4°. 

tion of gradients in WENSS LSS in Fig.|3lcan be identified as 
the two modes of the bimodal WENSS distribution shown in 
the bottom panel of Fig.|3] 

When ppfac increases, the BS signal becomes stronger 
compared to the WENSS signal, and the area where the BS 
signal is significantly stronger than the WENSS signal also in- 
creases. Both factors will 'convert' gradients from the steep 
gradients in the peaks in the bimodal WENSS distribution to 
flatter BS gradients that are found in the peak around radi- 
ans/deg. 

5.4. Dependence on interval length 

In Sect. l5.2l we discussed the distribution of the fitted gradients 
along Galactic latitude d(ppo\/db for an interval length of 1.4°. 
The distributions in the top panels of Fig.|3]are clearly bimodal, 
with a central dominant peak around radians/deg and a sec- 
ond clear peak at about -2.3 radians/deg. Furthermore there is 
tentative evidence for a plateau or small peak at around -1-2.5 
radians/deg. For the shortest 2 interval lengths we probed (0.7° 
and 1.1°) only the central peak is visible, with a plateau extend- 
ing to about -2.3°/bin. This plateau turns into a clear peak for 
intervals longer than 1.4°. At longer intervals (2.1°, 2.8°, 4.2° 
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Fig. 4. Average gradient per strip along Galactic latitude as a 
function of Galactic longitude. In the top panel we plotted gra- 
dients fitted to datapoints where the polarized intensity from 
WENSS is less than 3 times the polarized intensity in BS, in the 
bottom panel we plotted gradients fitted to the remaining dat- 
apoints. In the lower panel we show the average gradient per 
strip for the positive and negative modes separately. We only 
plotted strips that had at least 4 usable gradients, in the bottom 
panel each mode has at least 4 usable gradients. The errorbars 
are calculated from the spread (standard deviation) of the gra- 
dients in a strip, and are at the Icr level. The weighted average 
of the datapoints is indicated by a solid line, and in the bottom 
panel the weighted average was calculated for the positive and 
negative modes separately. The width of the interpolation ker- 
nel was 5.7° FWHM, ppimc-^, and the interval length of the 
fitted gradients is 1.4°. 



and 5.6°) the width of the central peak decreases and the pos- 
itive peak at +2.5 radians/deg becomes more prominent. For 
interval lengths of 4.2° and 5.6° the number of gradients in the 
central peak decreases rapidly. 

The decreasing width of the peak of gradients around 
radians/deg when going to longer intervals can be explained 
by the fact that longer intervals will span positive and neg- 
ative gradients on shorter lengthscales, and these 'subgradi- 
ents' will partially cancel each other At the same time the bi- 
modal WENSS distribution is found for all interval lengths in 
the WENSS-only data. These effects combined account for the 
positive WENSS mode becoming more pronounced when go- 
ing to longer intervals. 

The average x^^^ of the fitted gradients increases when we 
use longer intervals, going from < 4 for interval lengths of 0.7° 
to > 20 for intervals of 5.6°. This is due to structure in po- 
larization angle on scales smaller than the length of the fitted 
gradient. 
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Fig. 5. Average gradient per strip along Galactic longitude, as a 
function of Galactic latitude, plotted on the same vertical scale 
as in Fig. 12 We only considered strips that have a least 4 usable 
gradients. The errorbars are calculated from the spread (stan- 
dard deviation) of the gradients in a strip, and are at the Icr 
level. FWHM of the kernel used in the reconstruction is 5.7°, 
PPfrac=l, interval length over which the intervals are fitted is 
2.8°. 

6. Quantifying the large-scale stratification of 
polarization angle: results 

6.1. Gradients along Galactic latitude 

Fig. 0] shows the average polarization-angle gradients along 
Galactic latitude d(ppo\/db for each Galactic longitude in our 
sample, where we separately fitted gradients to datapoints 
where the ratio of polarized intensities Pwenss^Pbs is < 3 or 
> 3 (top and bottom panels respectively). This neatly separates 
the two distribtutions of slopes. In the bottom panel of Fig.l^we 
furthermore calculate the average gradient per strip separately 
for the positive and negative modes of gradients. Gradients in 
the top panel do not show this bimodal distribution. The solid 
line shows the weighted average over all Galactic longitudes; 
in the bottom panel the weighted average was calculated for the 
positive and negative modes separately. In both figures we only 
considered strips that have at least 4 usable gradients, a mini- 
mum number for an average gradient per strip to be statistically 
reasonable. 

From Fig. |2 it is clear that the gradients in the areas in 
WENSS LSS where WENSS is relatively bright or faint have 
very different distributions: whereas the WENSS faint areas 
show no clear gradients, the WENSS bright areas show a clear 
bimodal distribution of gradients at about +2 radians/deg and 
about -2 radians/deg over a wide range of Galactic longitudes. 
Furthermore this bimodal distribution is present for all the in- 
terval lengths we probed (from 0.7° to 5.6°). This last fact leads 
us to conclude that polarization-angle structure on scales that 
are small compared to the interval length do not dominate our 
gradient-fitting analysis (these structures only increase the 
of the fit), and that we are indeed sampling an underlying large- 
scale distribution. 



6.2. Gradients along Galactic longitude 

Gradients along Galactic longitude d(f)poi/dl behave different 
from gradients along Galactic latitude. In Fig. |5] we show 



Schnitzeler et al.: The WENSS & Dwingeloo surveys and the Galactic magnetic field 



9 



RS bright 




160 150 140 
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Fig. 6. RM and direction of Bj_ derived from BS polarization- 
angle measurements in the area covered by WENSS (enclosed 
by the dotted line) and its surroundings. The BS bright areas are 
shown as grey filled polygons. RM are shown when the fit to 
the observed polarization angles has;ifjgj < 5, and polarization 
angles should be available for at least 4 of the 5 wavelengths. 
The scale of the plotted symbols is linear and is illustrated in 
the top panel, units are radians/m^. 



the average gradient per strip along Galactic longitude. We 
only plotted strips that contain at least 4 usable gradients. The 
weighted average of the datapoints is indicated by a solid line. 
Since there is no bimodal distribution of the gradients we did 
not have to separate the fitted gradients as we did for gradi- 
ents along Galactic latitude. From this figure it is clear that the 
average gradient per strip is not significantly different from 
radians/deg. In the original WENSS dataset we found the same 
result for gradients along Galactic longitude. The increased 
spread of the average gradient per strip and also the increased 
errors per datapoint above a Galactic latitude of a: 18° coin- 
cides with a change in scale of the structure in polarization an- 
gle, which we noted earlier in the WENSS LSS mosaic in Fig. 
n and/or a decrease in signal-to-noise ratio. 

For the different interval lengths that we probed for strips 
along Galactic longitude (1.4°, 2.1°, 2.8°, 4.2° and 5.6°) the 
distribution of the gradients remains single-peaked, and the av- 
erage gradient along Galactic longitude is always close to 
radians/deg, but the average x^^^ per gradient increases as the 
interval length increases from < 5 to about 10 for the longest 
intervals. 



7. Comparison to previous results 

We assume that the amount of Faraday rotation is given by 
- 00 = RMA^. This means that if there are no gradients 
in the intrinsic polarization angle of the emitted radiation 
the gradients in polarization angle d(ppo\/dx (where 'jc' is either 
Galactic longitude or latitude) can be directly translated into 
gradients in RM with: 
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Fig. 7. Histograms of the gradients fitted to the RM we derived 
from the BS polarization-angle measurements (top row), and 
the gradients (l/A^) ■ dcp^oi/db derived from WENSS LSS (bot- 
tom two rows). The RM gradients we derive only use RM with 
^red ^ ^' ^^'^ polarization angles should be available for at least 
4 of the 5 wavelengths. The bottom row shows gradients in 
the WENSS LSS data after we convolved it with a 2.3° beam. 
Gradients in the left column were fitted to the BS faint data, 
gradients in the right column were fitted to the BS bright data. 
The gradients fitted to the BS RM were at most 5° long, the 
gradients fitted to the WENSS LSS data were 4.2° long. See 
the text for details. 

<po varies very slowly with position in the region we are consid- 
ering (see the directions of B± plotted in Fig.|6}, which means 
that the contribution of d(p{)ldx to dRM/dx is small. When we 
include the measured variations in 0o in our analysis we find 
identical results. 

To test the reliability of our 'translation' of polarization- 
angle gradients as RM gradients we calculated gradients in RM 
based on the BS dataset. We rederived the RM from the BS 
polarization-angle measurements to get also a x^^^ measure of 
the quality of the fitted RM, using the 1/signal-to-noise esti- 
mate of the error in the measured polarization angle described 



dRM 1 d(p. 



pol 



dx dx 



(4) 



in Brouw & Spoelstra (119761 . In calculating RM gradients we 
only considered RM with xl^.^ < 5, based on polarization an- 
gles for at least 4 of the 5 observing frequencies. In Fig.|6lwe 
show these RM. The RM gradients we calculated from the BS 
data use pairs of RM for which the separation in Galactic lon- 
gitude was smaller than 0.2 times the separation in Galactic 
latitude to constrain the orientation of the gradients fitted to the 
RM. We only fitted gradients to lines-of-sight that were less 
than 5° apart, which is slightly larger than the average distance 
between the datapoints shown in Fig.|6l 

Spoelstra ( 1984 1 points out that for a given sampling of 
space the smallest distance between two consecutive data- 
points will set a maximum to both the RM that can be deter- 
mined and to the B\\ that is derived from this RM. Spoelstra's 
argument is however based on only 2 datapoints, and the other 
available datapoints will make this criterion less strict. The RM 
that were calculated by Haverkorn et al. (i2003aj are compara- 
ble in size to the RM found by BS in the same region, but are 
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derived from 5 MHz wide bands around about 350 MHz. An 
additional 180° difference in polarization angles between two 
such consecutive bands would require a RM that is about 150 
radians/m^ larger than they derive. In the direction we are con- 
sidering such large RM are only observed in pulsars that lie in 
or beyond the Perseus arm, and since depolarization effects will 
limit the line-of-sight, we do not believe that such large RM 
should be observable in these data. Therefore we think that also 
the BS RM are coiTect and do not suffer from 'aliasing' due to 
the poor sampling. 

In Fig. we show the gradients we fitted to the BS RM 
separately for the BS faint and bright areas. If we found a sig- 
nificant difference between the polarization-angle distribution 
in WENSS LSS and the original distribution in WENSS such 
an area is identified as BS bright. In Fig.|6lwe indicated the 
BS bright areas by a grey shading. In the second row of Fig. 
we show the distribution of the gradients we fitted to the 
polarization-angle data, scaled to gradients in RM by the factor 
of 1 lA^ in Eqn. |2 These gradients were fitted over an interval 
of 4.2°. 

To see if the distribution of the BS RM gradients can be 
scaled up to the gradients we derived from WENSS LSS we 
used a Kolmogorov-Smirnov two sample test using Eqn. |3 as 
a description for the mapping of the BS gradients onto the 
WENSS gradients, where we apply a scale factor to the BS 
RM. At the 5% confidence level for rejection scale factors be- 
tween 1.7 and about 3.1 are permitted, with an average value 
of about 2.4. The KS probability has a maximum of 0.7 for a 
scale factor of 2.2. 

In the BS faint region the RM gradients we derived are 
about a factor of 2.4 smaller than the RM gradients implied 
by the WENSS data. To simulate the effects of the difference 
in size between the WENSS and BS beams we convolved the 
WENSS LSS map with the 2.3° beam of the BS data at 408 
MHz. The results for the convolved data are shown in the bot- 
tom row of Fig.0 Clearly the larger telescope beam of the BS 
data can account for the narrowing of the distribution and the 
apparent disappearance of the negative mode. 

8. Physical picture 

In this section we will interpret our data in terms of a simple 
model of the ISM. In Sect. [O] and lO we will derive infor- 
mation on the magnetic field component along the line-of-sight 
and on the component perpendicular to the Galactic plane. In 
Sect. l8.3l we will discuss what we infer from our data about the 
so-called 'fan' region. 

In our model we simulate a single line-of-sight through a 
volume of ISM with synchrotron emission and/or Faraday ro- 
tation. This volume contains both a regular and a random mag- 
netic field. The regular field has a constant strength and direc- 
tion, and the random field has a constant strength, but it changes 
direction on a scale that is (much) smaller than the size of the 
simulated volume. This ensures that there are enough draws of 
the orientation of the random field component along the line- 
of-sight that its statistical properties can be calculated analyti- 
cally. The amount of emission only depends on the length of the 
total magnetic field vector (the sum of the regular and random 



components) projected perpendicular to the line-of-sight, as in- 
dicated in Appendix A of Haverkorn et al. (2004i. Similarly, 
the total magnetic field projected along the line-of-sight sets 
the amount of Faraday rotation. However, only a fraction of the 
line-of-sight is filled with cells that contain Faraday rotating 
electrons; in each of these cells the electron density is constant. 
The simulated volume of ISM is illuminated from the back by 
polarized synchrotron radiation; this background is assumed to 
have no internal Faraday rotation. 

One important ingredient in this model is the dispersion 
measure DM, the line-of-sight integral of the electron den- 
sity, DM = / n,d/, [DM] -cm ^pc, which we cannot deter- 

l.o.s. 

mine from our observations. To estimate DM we integrated the 
NE2001 electron density model by Cordes and Lazio (,2002> 
out to 500 pc ('local') and 2 kpc (about the distance of the 
Perseus arm for these lines-of-sight). Since pulsars in the sec- 
ond Galactic quadrant that lie beyond the Perseus arm have RM 
that are typically 100 radians/m^, at least an order of magnitude 
larger than RM of pulsars that lie in front of the arm (see the 
Galactic distribution of pulsars in e.g. Weisberg et al. 200?t . 
we argue that the magnitude of the BS RM, typically less than 
10 radians/m^, puts the Faraday -rotating electrons that produce 
these RM between us and the Perseus arm, which means that 
the 2 kpc we use is a reasonable upper limit to the length of the 
line-of-sight. Note that for the 2 kpc line-of-sight and a Galactic 
latitude of 18° (the maximum latitude that we will probe - see 
Table|2}, the line-of-sight reaches a maximum height above the 
Galactic plane of about 600 pc, which means that we do not 
reach out into the Galactic halo. 

8.1. Deriving the iine-of-sight magnetic field 
component 

In this section we will derive strengths for the magnetic field 
component along the line-of-sight for a number of Galactic lat- 
itudes. We will use the RM determined from the BS data and 
the DM we calculated from the NE2001 model for the 2 as- 
sumed lengths of the line-of-sight. 

Following Sokoloff et al. (1998 1 we introduce the 'intrin- 
sic rotation measure' fi defined as the line-of-sight integral 
of the product of the magnetic field and the electron den- 
sity = 0.81 / rieB ■ dl, [B]=juG, [«,]=cm-3, [|/|]=pc, and 

l.o.s. 

["^^l^radians/m^. As is customary the length of a vector will be 
denoted by Z? = \B\. Since it is an integral along the line-of- 
sight it will not depend on the 'dumpiness' of the distribution 
of the Faraday rotating electrons. 

If the line-of-sight is long enough, or if the average is taken 
over a number of lines-of-sight that are short, the contribution 
of the random magnetic field to will cancel out in our model, 
and 'R oc DMZ?reg,||- This means that we can calculate fireg.n from 
H and DM using 

observer 

J HeB ■ dl 

= QMDM = nAl " 

l.o.s. 
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where the magnetic field strength is in fiG. The ratio of integrals 
on the r.h.s. of Eqn. |5] defines the electron density weigthed 
average -Bi-eg,||, or (-Sreg,||)n,- Since the -Breg.y we calculate from 
Eqn.|5]is defined to be constant throughout our model, (Breg,||)«<, 
will be equal to Breg.n in this case. 

The rotation measure RM is defined as the derivitive of po- 
larization angle with respect to wavelength squared and, con- 
trary to H, it does depend on the exact distribution of emitting 
and Faraday rotating regions along the line-of-sight. This has 
been shown for example in Sokoloff et al. ( 1998 1. Therefore the 
ratio of RM to will be different for different configurations 
of emitting and Faraday rotating regions. 

We use BS polarization-angle measurements with Galactic 
longitudes between about 150° and 160° and Galactic latitudes 
between 5° and 20°, excluding the Auriga region between 156° 
< / < 165° and 12° < < 18° where the RM might be affected 
by local structure. Most of the BS data we consider thus lie in 
an area of WENSS LSS where LSS does not dominate. We es- 
timated RM and the intrinsic polarization angle of the emitted 
radiation 0o from neighbouring (measured) lines-of-sight along 
6 Galactic latitudes at a Galactic longitude of 152°, as shown in 
the first 3 columns of Table|2l and the errors given in the head- 
ings of these columns reflect our estimate of the spread around 
these averages. 0o can be derived by extrapolating the BS po- 
larization angles to A-Q meters. Note that in the 3rd column 
of Table |2] we indicate the direction of the magnetic field per- 
pendicular to the line-of-sight, B^, which is perpendicular to 
the direction indicated by ^q. At other Galactic longitudes the 
number of reliable RM is either much lower, or LSS is miss- 
ing from the interferometer observations. To correct the RM in 
Table|5]for beam effects these values had to be multiplied by a 
factor of 2.4 + 0.7 as discussed in Sect.0 

In columns 4 and 6 of Table |2l we show the values for 
^reg.ii we calculated from Eqn. [S] for the two assumed lengths 
of the line-of-sight. The error in Bieg.n is based on the range 
in Galactic latitude of the position of the line-of-sight (which 
gives a range in allowed DM), the error in the estimated RM, 
and the uncertainty in the scale factor between the WENSS and 
BS RM gradients. We compared the Breg.n found using Eqn.ls] 
to the Breg.ii from the model by Haverkorn et al. J2004> . in the 
remainder referred to as 'H04'. In their model Faraday rota- 
tion occurs in the synchrotron-emitting thin disc in cells filled 
with thermal electrons. Cell sizes range between 1 and 60 pc. 
These cells have a constant electron density of 0.08 cm"^, and 
the filling factor assumed in this model is 0.2 (values based 
on Reynolds I1991I I. The height of the synchrotron emitting 
thin disc is 180 pc (Beuermann et al. I1985l l. The H04 model 
contains both a regular magnetic field, with constant strength 
and direction, and a random field, with constant strength, but 
a direction that is randomly drawn in each cell. The thin disc 
is illuminated from the back by a synchrotron-emitting halo. 
Emission and Faraday rotation in this halo are assumed to oc- 
cur on such large angular scales that the emission from the halo 
has a uniform total and polarized intensity and direction of the 
plane of polarization. 

The H04 model simulates a large number of lines-of-sight 
that each have their own distribution of Faraday rotating cells. 
By specifying the strength and the direction of Breg.n 



Breg,±, the strength of Bran, and the total and polarized intensi- 
ties of the background illuminating the thin disc, the H04 model 
predicts by numerical radiative transfer the observables Stokes 
1, RM, Stokes Q and U and their standard deviations (the spread 
of the observables being created by the different configurations 
of Faraday rotating cells along the line-of-sight). By comparing 
these output parameters to their observations, H04 constrain the 
range of allowed input parameters. 

One major advantage of the H04 model over our own ana- 
lytical model is that the numerical treatment in the H04 model 
can predict RM. RM is difficult to calculate analytically, since 
every configuration of Faraday rotating cells will have a differ- 
ent amount of depolarization, which means that the behaviour 
of the polarization angles with wavelength squared will be dif- 
ferent. Therefore our model depends on a relation between RM 
and as input, which we inferred by comparing our results to 
the ones from the H04 model. 

H04 find a Breg.n of -0.42 yuG in the direction of their Auriga 
field. When we use Eqn.|5lwe find a value of -0.40 fiG if we 
use the same input DM as in the H04 model, and by setting 
K - RM. From this we conclude that setting K - RM works 
well in this case. The DM used in the H04 model in this direc- 
tion is 10 cm^^pc, lower than the 3 1 cm^^pc we determine from 
the NE2001 model by integrating it out to the same distance. In 
part this difference in DM is due the filling factor that is used in 
the H04 model. For the clumpy disc they model, where all cells 
have the same electron density, a filling factor of 0.4 instead 
of the 0.2 used in the H04 model would agree better with the 
Reynolds data. This doubles the DM in the H04 model, bring- 
ing it closer to the value predicted from the NE2001 model. 

8.2. Reconstructing the 3D magnetic field vector. Tlie 
magnetic field strength perpendicular to the 
Galactic plane 

In this section we will derive information on the magnetic field 
component perpendicular to the Galactic plane B^. Since our 
derivation will only apply to the regular magnetic field we will 
drop the 'reg' subscript when we indicate field components. 
The relations between the different magnetic field components 
are illustrated in the top panel of Fig. |S] First we will derive 
the strength of the magnetic field component perpendicular to 
the line-of-sight B± and its orientation. Bj_ can be decom- 
posed into a component along Galactic longitude and a com- 
ponent along Galactic latitude By,. The component of Bj_ along 
Galactic longitude does not contribute to B^. B^ can then be 
calculated from the projected B^ and B\\: - B^^^ + B~^\\ - 
Bb cos(b) - B|| sin(|/7|), as illustrated in the bottom panel in Fig. 
|8] The sign of B. is such that a positive B~ points away from the 
Galactic plane both for positive and negative Galactic latitudes. 
The sign of Bi, should be chosen accordingly. 

The strength of Bj_ cannot be determined from our data. 
However, if we interpret the different By that we determine in 
Table|2las the result of different orientations of B with respect 
to the line-of-sight (assuming that B has a constant length for 
these lines-of-sight), B^ can be calculated as a projection of B. 
To do this B has to be known. In the vicinity of the sun pulsar 
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Observer 



sideview (in a plane of constant Galactic longitude} 

Fig. 8. Scheme relating the different magnetic field components 
that are discussed in the text. Panel B shows panel A from the 
side in a plane of constant Galactic longitude. The magnetic 
field components along and perpendicular to the line-of-sight 
are B\\, resp. B^. The plane spanned by B\\ and Bj. is indi- 
cated in grey. Bj. is directed perpendicular to both being 
defined in the plane perpendicular to the line-of-sight. Bj, is 
the component of B^ along Galactic latitude. We indicated the 
directions of increasing Galactic longitude / and latitude b and 
height above the Galactic plane z, and defined Bt, and B, to 
increase when moving away from the Galactic plane. 



RM give a B that is typically about 2 fiG, and equipartition ar- 
guments indicate for the total magnetic field (regular +random) 
a strength of about 6.5 //G (Beck et al. ri996t . In H04 both By 
and Bj^ are derived, which give a B ^ 3.3 juG. If we assume 
that B in the direction of H04's Auriga field is the same as in 
our data, we only have to rescale their B to the lengths of the 
lines-of-sight we use. The scaling of B|| depends only on DM, 
and we can use our analytical model to derive the scaling prop- 
erties of Bx from Stokes I and ctj, in a similar way as was used 
in H04. In this way we find B = 3.2 + 0.5 and 2.0 ± 0.4 //G for 
our 500 pc and 2 kpc lines-of-sight. 

B/, can be derived from B/, - B^ sin(2(0o -i- 7t/2)). Note that 
since (po follows from 1/2 arctan(U/Q), where U and Q refer 
to Stokes U and Q resp., it has a period of 180°, which is why 
there is a factor of 2 in this equation. The plane of polarization 
(and therefore also the direction of B±) is perpendicular to ^o, 
which accounts for the n/l in the expression for Bh- Once the 
orientation of the plane of polarization is known, Bj_ can still 
point in 2 directions. The global geometry of the magnetic field 
can be used to determine in which of the two possible direc- 



tions Bj_ is pointing: Brown et al. ( 12003> and Johnston-HoUitt 
et al. J2004t derive from RM of extragalactic radiosources that 
the field points towards smaller Galactic longitudes. We will 
assume that the features we observe in WENSS LSS follow the 
same global field direction. 

The errors we quote for By and B, are based on the range in 
Galactic latitude for the position of the line-of-sight, the spread 
in RM and the orientation of and (in the case of B,) also 
uncertainties in the parameters modeled by H04. Even though 
our estimates of B; are Icr values at best, Br is between 1/2 to 
4 times By, for both lengths of the line-of-sight that we use. 
Since in all these cases the By and B/, components are small 
compared to the total B, most of B should be oriented along 
Galactic longitude, and therefore the magnetic field must be 
azimuthal, and directed towards smaller Galactic longitudes. 
Using pulsar observations Han and Qiao ( .1994J estimate that 
B, ^ 0.2 - 0.3 yuG, similar to our results. 

To reach these conclusions we made a number of assump- 
tions that, though reasonable, are hard to prove using our cur- 
rent data. We did not take into account the intrinsic difficulty of 
constraining the NE2001 model in the second Galactic quad- 
rant which is the result of the small number of known pulsars in 
this direction. Since this only influences DM and therefore By, 
but does not influence Bj^, we think it unlikely that the error in 
DM is so large that it influences our conclusion of an azimuthal 
field. With better datasets, of the type described in Haverkorn 
et al. (2003aland l2003bl l. some of these assumptions could be 
tested for their validity, and the errors in By and Bj. could also 
be reduced. We have obtained several of such datasets recently. 

To investigate if the results we present here can be consid- 
ered 'typical' for this quadrant we can look at the BS RM that 
lie outside the area covered by WENSS, since we showed in 
Sect, □that the BS RM can be considered as scaled-down ver- 
sions of WENSS RM (at least in the BS faint part of WENSS 
LSS). Outside WENSS the BS sampling however becomes 
much poorer, and also the RM seem to become much less 
correlated (which is probably in part due to the poorer sam- 
pling). Therefore it is at the moment not possible to investigate 
whether the results we present here apply only to a local feature 
or are more general. 

8.3. Implications for the fan region 

The region between 136° < I < 152° and -3° < b < 13° is 
part of the so-called 'fan' region that we already mentioned in 
Sect.0] In the original WENSS dataset (without the single-dish 
information) the polarization-angle gradients in this area are bi- 
modally distributed like in the rest of WENSS, as is also clear 
from the bottom panel in Fig.|3 where there is a lack of gradi- 
ents around radians/m^/deg even though this figure includes 
polarization-angle gradients in the fan region. One difference 
between the fan and the region outside the fan is that the ratio 
of positive to negative polarization-angle gradients in WENSS 
inside the fan is about 1:1, whereas this ratio is more 1:3 for 
WENSS as a whole. 

We suggest that the fan is a region of enhanced emission 
that lies in front of the emission we see in WENSS. Due to 
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Table 2. Estimated average RM and direction of Bj^ ('Bj^') from the BS data at a longitude of about 152° for different Galactic 
latitudes. Estimates of the uncertainties in the parameters are given in the fourth row and in the columns for Bieg,||, and are based 
on the range in Galactic latitude of the position of the line-of-sight, the spread of both RM and direction of B^, the uncertainty 
in the scale factor between the WENSS and BS RM gradients, and the modeling uncertainties in the H04 model. Using Eqn.|5] 
the measured RM (after being corrected for beam effects as described in Sect.^} together with DM calculated from the NE2001 
model out to 500 pc and 2 kpc, can be used to calculate the magnetic field component along the line-of-sight, Breg.n- B. can then 
be calculated by projecting B\\ and the component of Bj_ along Galactic latitude, as is illustrated in the bottom panel of Fig.lH 
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* the b = 1° datapoint lies in a part where the BS data is 'intermediate' in brightness, whereas all other datapoints are in BS faint areas. The 

B could therefore be different compared to other datapoints. 
t at these latitudes 0o could not be reliably estimated, which means that could not be determined. 



the short lines-of-sight to the fan in this model, Faraday rota- 
tion effects, in particular variations in the amount of Faraday 
rotation across the fan surface, will be small. The polarized 
signal emitted by the fan will then be modulated by the fore- 
ground Faraday screen on angular scales that are too large 
to be detected by an interferometer due to its missing short 
spacings. The longer lines-of-sight to the background emission 
would make Faraday rotation effects strong enough that the 
background can be detected by the interferometer We there- 
fore identify structure we see in WENSS as resulting from 
the Faraday modulation of this background emission. If the 
fan would lie at the same distance as the WENSS emission it 
should also be visible in the interferometer data. In the large 
single-dish beam the background emission, which is highly 
structured in polarization angle due to Faraday rotation, will 
get beam depolarized. If the fan would lie at the same distance 
as the background it too would get beam depolarized in the 
single-dish data, contrary to what we observe. 

In this model the separation along the line-of-sight be- 
tween the fan and the background emission is determined by 
the amount of small-scale structure in polarization angle. If 
the amount of structure on small angular scales is smaller than 
what we derive from WENSS (e.g. due to missing short base- 
lines) the separation in distance becomes less. However the 
similarity of BS RM inside and outside the fan, supplemented 
by the presence of a bimodal WENSS distributed in both these 
regions, suggests that we are not grossly overestimating the 
variation of polarization angles on small scales. 

Spoelstra (il984i) concludes from a comparison between the 
polarization angles of starlight emitted by stars at known dis- 
tances and 00 that the Faraday rotating medium responsible for 
the RM observed towards the fan should lie within a few hun- 
dred parsecs. RM derived from the diffuse ISM in the direction 



of the fan are small compared to RM derived for pulsars that, 
in this part of the Galaxy, lie beyond the Perseus arm, which 
would also indicate that most of the Faraday rotating medium 
seen towards the fan lies in front of the Perseus arm. Wolleben 
(2005) suggests that the source of the emission seen towards 
the fan region forms a single physical structure that extends 
both above and below the Galactic plane. He tentatively iden- 
tifies the low polarized intensities near the Galactic plane that 
appear to separate the parts of the fan region above and be- 
low the Galactic plane as depolarization occuring in HII com- 
plexes. This puts the source of the emission observed towards 
the fan beyond these complexes. One of the furthest of these 
complexes (with a moderate coiTelation with polarized inten- 
sity) is the IC1795/1805 complex around {I, /7)=(135M°) at 2. 1 
kpc. The observations by Wolleben et al. were caiTied out at 1 .4 
GHz, whereas Spoelstra uses also data at much lower frequen- 
cies. Since depolarization (both line-of-sight and beam) wiU be 
more important at lower frequencies, the bulk of the observed 
signal will be coming from closer to the observer than at higher 
frequencies, which makes comparing data at such different fre- 
quencies more difficult. 

The fan region can clearly be seen in the RM in Fig. |6l 
One striking feature appears to be a (smooth) change in RM 
from positive at around {I, b)-{l50°,6°), to about at Z=142° to 
negative at around {l,b)-{\'iQ° ,\\°). Furthermore the negative 
and positive RM are of about the same magnitude. A physical 
model of this structure would require knowledge of the strength 
of the magnetic field component in the plane of the sky. The ob- 
servations and modelling by Haverkorn et al. ( »2003bJ imply a 
strong (si 3.2 juG) field in the plane of the sky, but their mod- 
elling includes also the 'ring' feature at {l,b)-(\?>l° ,1°). There 
are a number of indications that this ring structure could be a 
different type of object than its surroundings : the single large 
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RM in Fig. |5] lies at the position of the ring, whereas the other 
RM in its vicinity are smaller, and furthermore the single-dish 
data show a local depression at the position of the ring in the 
generally high polarized intensities seen towards the fan region. 

9. Summary & Conclusions 

We have studied the polarization properties of the diffuse 
Galactic radio background at 325 MHz in the region 130° < 
/ < 173°, -5° <b<3r using the WENSS dataset. We de- 
termined gradients in polarization angle (which can be trans- 
lated into gradients in RM), taking into account the periodic 
nature of the polarization-angle data. The WENSS data show 
gradients that are coherent over large areas. We find a bimodal 
distribution of gradients along Galactic latitude, with peaks at 
+2.1 radians/deg and -2.1 radians/deg, with about three times 
as many gradients in the negative mode as there are in the pos- 
itive mode, and with hardly any gradients in between the two 
modes. Along Galactic longitude the gradients are in general 
nearly zero. 

We investigated the importance of LSS missing in our 
data from the higher-frequency single-dish data obtained by 
Brouw&Spoelstra (I1976t . We first interpolated between the 
lines-of-sight of the single-dish data, followed by an extrap- 
olation to 325 MHz, the observing frequency of WENSS. We 
ran several consistency tests between the different frequencies 
present in the single-dish observations which show that the data 
allow this handling in a robust way, independent of the exact 
parameters we use. Our reconstruction of the LSS missing from 
WENSS is only statistical in nature, we do not claim to recon- 
struct the correct observables for individual pixels. We refer to 
the combined WENSS/BS datasets as 'WENSS LSS'. We es- 
timate that the combined dataset is still insensitive to structure 
on spatial scales between ail. 5° and a;5°. 

After adding our estimate of missing LSS to the WENSS 
data, many polarization-angle gradients are 'converted' from 
the steep bimodal WENSS distribution into a peak around zero 
radians/deg. The negative mode of WENSS gradients remains 
clearly visible. Since we constructed WENSS LSS as the (po- 
larization) vector sum of WENSS and BS, the original WENSS 
or BS data will dominate if its polarization vectors are much 
longer than the vectors in the other dataset. Therefore in re- 
gions where the BS signal is strong the slowly-varying (nearly 
zero radians/deg) BS gradients will dominate WENSS LSS. 

By using a Kolmogorov-Smirnov (KS) test to compare the 
polarization-angle gradients from the BS faint part of WENSS 
LSS to scaled-up gradients from the BS RM. We found that 
these distributions are similar, and we set limits to the scale 
factor using this KS test. We illustrated the importance of 
the difference in size between the WENSS and BS beams 
by convolving the WENSS LSS data with the BS 408 MHz 
single-dish beam, which made the clear bimodal distribution 
of polarization-angle gradients in WENSS LSS look more like 
the distribution of BS RM gradients. 

We proceeded by deriving the strength of the magnetic field 
component along the line-of-sight using these scaled-up RM 
together with dispersion measures from the NE2001 thermal 
electron density model by Cordes & Lazio (I2002> . see Table 121 



where we show the values derived for lines-of-sight of lengths 
0.5 and 2 kpc. By comparing the magnitude of the scaled-up 
RM to RM measured for pulsars in front of and beyond the 
Perseus arm we conclude that the Faraday rotation effects we 
are observing in the BS faint part of our data should be occuring 
between us and the Perseus arm. 

We used the estimate of the strength of the total (regular) 
magnetic field by Haverkorn et al. (2004i to derive the strength 
of the magnetic field in the plane of the sky Bj_, and even 
though our results are about l-3cr at best, we showed that the 
component of B± along Galactic longitude is much stronger 
than both the component of B± along Galactic latitude and B\\. 
We were also able to determine the strength of the magnetic 
field component perpendicular to the Galactic plane, B., shown 
in Table 12 

Part of the area covered by WENSS LSS includes the 'fan' 
region, an area bright in polarized intensity, but at the same 
time extended on large angular scales, which meant that it 
could not be detected by the WSRT. We note that the large- 
scale emission that is detected by the single-dish does not vary 
on the same angular scales as the WENSS data, and we sug- 
gest that the large-scale emission seen towards the fan lies in 
front of the emission we see in WENSS. The magnitude of the 
RM observed towards the fan are small, indicating that they 
are caused by Faraday rotation occurring between us and the 
Perseus arm. 
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Appendix A: A new method for fitting linear 
gradients to periodic data 

The main problem when fitting functions to periodic data like 
polarization angles is the correct handling of the 180° ambi- 
guity. Finding the best-fitting slope and offset requires at the 
same time finding the number of times 180° has to be added to 
each datapoint to minimize the distance between the Une fitted 
to the data and the datapoint. This means that for n datapoints 
n + 2 degrees of freedom have to be solved, 1 degree of free- 
dom per datapoint to determine the number of 180° 'wraps' 
(A^,; i=l,...,n), and 2 degrees to determine slope and ofi'set. The 
only solution to this problem is that at least 2 of the A',- are 
known a priori. The best-fitting slope and offset will however 
only depend on the (small number of) datapoints for which the 
A', are known, which can be shown by working out the equa- 
tions that minimize x^- 

An alternative approach which takes the behaviour of all 
datapoints into account when determining slope and offset can 
be formulated as follows. When making a linear fit one can pick 
a combination of slope and offset and then bring each datapoint 
to within 90° of the imposed linear fit, which will minimize the 

for that fit. By comparing the x^ for all possible combina- 
tions of slope and offset one then finds the best fit to the data. 

Instead of considering all possible slopes and offsets gener- 
ally only combinations are considered that he on a finite grid. 
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Fig. A.l. Illustration of which domains and configurations are 
allowed by a simple example distribution of 4 pixels with co- 
ordinates (pixel,polarization angle)= (1,0°), (2,45°), (3,-90°) 
and (4,-45°). All domain boundaries that belong to a single 
datapoint run parallel to each other, and the boundaries are 
steeper if the datapoint is further from the origin of the fitted 
line, which is at pixel coordinate in this example. As a com- 
parison we also show a grid that probes (slope.offset) combina- 
tions (indicated as crosses) with grid spacing 20°/pixelx20°. 



For a periodic function only the central part of the complete 
(slope,offset) space has to be investigated. For polarization- 
angle gradients this means that only slopes between -90°/pixel 
and +90°/pixel and offsets between -90° and H-90° have to be 
investigated since other slopes and offsets are indistinguishable 
from slopes and offsets in this central space. 

In this appendix we discuss an alternative method which is 
not based on a grid of slopes and offsets. Instead it maps areas 
in (slope,offset) space ('domains') which share the same con- 
figurations (Ni,..., Nn) of 180° wraps. For a given (slope,offset) 
combination the lowest x^ is found when all datapoints are 
within 90° of the fitted line formed by that (slope,offset) combi- 
nation. Each domain boundary is thus formed by (slope,offset) 
combinations for which one datapoint is at exactly 90° from the 
fitted line. Crossing this boundary means going to (slope,offset) 
combinations for which the datapoint will be further than 90° 
from the fitted line, and it can be brought to within 90° by 
adding an extra + 180° to the datapoint, which changes the con- 
figuration. Each datapoint thus defines a set of domain bound- 
aries, and the ensemble of the domain boundaries of the indi- 
vidual datapoints will define the domains. 

For each of these domains the best-fitting combination of 
slope and offset is found by a x^ fit- Since all combinations 
of slopes and offsets in the domain share the same best fit, it 
becomes possible to cover the entire solution space by mapping 
the domains. 

One possible cause of concern could be that the least- 
squares solution falls outside its domain. However the x^ of 
this solution will always be higher than the minimum of the 
best fit in the domain the solution has crossed into: If the best- 
fitting (slope.ofFset) combination Ues outside the domain it is 
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supposed to be in, by definition of the domain borders at least 
one of the datapoints is further than 90° from the line formed 
by that (slope,offset) combination. For that (slope,offset) com- 
bination it is however always possible to find a configuration 
where all datapoints are within 90° of the line formed by that 
(slope,offset) combination, and this configuration will therefore 
give a better fit. This furthermore guarantees that the lowest 
in (slope,offset) space will always lie inside the domain it is 
supposed to be, and that by mapping the domains and using a 

fit in every domain we will find this best-fitting solution. 

The boundaries of the domains are set by the noise realiza- 
tion of the observations. Since the same distance between the 
datapoint and the linear fit is used to find the position of the do- 
main borders and in the definition of;^'^, the domain mapping 
approach will be no more sensitive to the influence of noise 
than a standard fit applied to a non-periodic dataset. 

The configuration-mapping approach covers all solutions 
in (slope,offset) space and is therefore much more reliable 
than the grid-based approach which only probes the quality of 
certain (slope,offset) combinations. But furthermore it is very 
easy to figure out where the boundaries of the domains are, 
and therefore this method is also much faster. In one example 
of 4 datapoints at (pixel,polarization angle)= (1,0°), (2,45°), 
(3,-90°) and (4,-45°) the grid-based approach using a coarse 
grid of 20°/pixelx20° already requires 100 gridpoints, whereas 
there were only 21 domains in this configuration that needed to 
be investigated (see Fig. lA.U . 

Sarala and Jain (2001 1 discuss a different method to find 
the best-fitting gradient for a periodic dataset. Their analysis 
includes all wraps of the individual datapoints and using the ap- 
propriate statistics they derive a maximum-likelihood criterion. 
By bringing the datapoints to within 90° of the fitted gradient 
our method can use the standard non-periodic statistics. 

The method we describe here can also be used for fit- 
ting RM to polarization-angle observations at different wave- 
lengths, since this would mean solving = aA^ + b for the 
parameters 'a' and 'b', which is identical to the (p - ax + h we 
solved for the spatial RM gradients. One important difference 
however is that the regular spacing of the pixels we fitted RM 
gradients to limits the range of slopes that give unique fits to the 
data: slopes steeper than 90°/pixel are indiscernible from slopes 
that are flatter If the polarization angles are irregularly sampled 
in space there are no such limits, which means that there are 
no boundaries that limit the range of RM that have to be in- 
vestigated. The regular spacing of datapoints also causes the 
degeneracy that steeper gradients can be fitted with the same 
quality by adding tt radians to the second datapoint. In to the 
third etc. For an irregular spacing this degeneracy is (at least 
partially) lifted. 
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